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Motivated by recent inelastic neutron scattering experiments on cuprates, we discuss the formation 
of bond order in the stripe phase. We suggest that the spin Peierls order emerges in hole-rich 
domain walls (DWs) formed between hole-poor regions in which long-range antiferromagnetic (AF) 
correlations exist. On the example of a single stripe we analyze the stability of such structures. 
The motion of a hole inside the DW which takes the form of a bond ordered ladder is in principle 
unrestricted. The hole hopping in domains is to some extent obscured by the fact that a moving 
hole spoils AF correlations. The propagation of a hole along the stripe which takes the form of 
the ladder-like domain wall that separates antiphase AF domains is a combination of these two 
types of motion occurring in two different environments. By analyzing the energy dispersion of 
a quasiparticle propagating along the bond-centered, bond-ordered stripe and of a quasiparticle 
propagating along the site-centered stripe we deduce that bond ordered stripes are stable at the 
total doping level 1/8 and the linear stripe-filling level 1/2. This conclusion seems to be relevant to 
the nature of the stripe phase in La\.s7z,Bao.\2z,Cu04,. 



PACS numbers: 74.20.Mn, 71.10.Fd 



I. INTRODUCTION 

The phenomenon of high temperature superconduc- 
tivity (HTSC) in quasi two dimensional copper oxides 
and other aspects of unconventional behavior observed 
in these systems are far from being understood by the 
community of physicists. Nevertheless much evidence has 
been got that the magnetism of these compounds influ- 
ences their transport properties. For example, an anoma- 
lous suppression of the critical temperature in supercon- 
ducting La2-xBaxCuOi and Lai,e-xNdoASrxCu04 at 
the doping parameter x in the vicinity of the commensu- 
rate level 1/8 1] has been associated with the formation 
of charge and spin order ^] . According to such a scenario 
which is based on neutron-scattering studies this ordering 
takes the form of stripes which may be viewed as charge 
filled DWs between AF regions that are formed in each 
copper-oxygen plane. That plane constitutes the most 
important building block in the crystallographic struc- 
ture of cuprates. The phase of the staggered magneti- 
zation in domains changes by tt across each DW. No di- 
rect measurements, which could determine the magnetic 
structure of DWs, have been performed. 

It is believed that the essential part of cuprate physics 
takes place in the single cooper-oxygen plane. Thus, we 
will concentrate on this structure from now on. There 
exists much evidence both experimental and theoretical, 
that the t-J model on the square lattice is well suited to 
describe the physics of the single copper oxygen plane. 
Axes of horizontal or vertical stripes which are regions 
with enhanced content of holes and which simultaneously 
form DWs may lead either through sites (in the case of so 
called site-centered stripes) or through centers of bonds 
(in the case of so called bond-centered stripes). In the 
framework of the theoretical model these sites and bonds 
form the square lattice at which the tJM is defined. 



It has been pointed out in recent papers Q that some 
high-energy excitations observed in neutron-scattering 
studies of Lai,g,75BaQ,i25Cu04 and YBa2Cu30e,e are 
distinct from conventional spin waves which propagate 
in the AF and are similar to excitations which are mea- 
sured in two-leg spin ladders. Much evidence has been 
obtained that the formation of the spin-Peierls order de- 
termines the physics of excitations in undoped and doped 
AF two-leg ladders P|. For example, an excitation which 
propagates in undoped system may be viewed as a bond 
triplet that moves in the background of bond singlets. 
The relevance of bond order, which takes the form of 
the spin-Peierls state, to behavior of doped 2D AFs and 
cuprates in particular, was suggested already some time 
ago 1^ . A class of states which have been discussed 



in this context comprises those having (Sj) = 0, because 
each spin finds a partner to form a bond singlet. Since 
there are many choices for the partner with which a given 
spin may pair up, the system may fluctuate between dif- 
ferent configurations. The theory of strong fluctuations 
between many conflgurations with different patterns of 
bond singlets covering the whole lattice, is called the res- 
onating valence bond picture 7]. In this paper we con- 
centrate on a scenario according to which fluctuations 
between different configurations are limited to small re- 
gions and a given pattern of bond ordering dominates. 
Theories of insulating states in cuprates with AF order 
in domains were being developed already in late 80-ties 
^j. Here we will try to answer the question if a uni- 
fying picture of coexisting magnetically ordered domains 
and DWs with the dominating spin-Peierls order provides 
the correct description of the stripe state in cuprates. 
In particular, we will try to check if the bond-centered, 
bond-ordered, hole-filled stripe is more stable than the 
narrow site-centered stripe. Later we will also analyze 
how doping and the formation of a stripe superlattice in- 
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fluences the outcome of the competition between these 
two types of structure. Through the paper we assume 
that the bond-centered stripe may be viewed as a DW 
which resembles a two-leg ladder filled with holes. 

The detailed plan of our analysis looks as follows. 
First we will discuss the motion of holes along the bond- 
centered ladder-like stripe with the spin-Peierls order ei- 
ther on rungs or on legs. Next we will calculate the energy 
of the "narrow" site-centered DW in the undoped quan- 
tum 2D AF and for the same system the energy of the 
"wide" DW which resembles the bond-ordered two-leg 
ladder. Finally, by calculating the energy at given lin- 
ear hole-filling we will compare the stability of the site- 
centered stripe with the stability of the bond-centered 
stripe in the 2D hole-doped quantum AF. We will also 
discuss how doping and the formation of a stripe array 
influences the underlying magnetic structure of stripes 
and their shape. 



II. COMPOSITE MOTION OF A HOLE ALONG 
A LADDER-LIKE STRIPE WITH SPIN-PEIERLS 
ORDER ON RUNGS 

We consider stripes which are formed in the doped AF 
between anti-phase domains as DWs filled with holes. 
It is not evident, at first glance, what is the magnetic 
structure of these DWs. The most obvious possible form 
of a single filled DW, which has been discussed in the 
literature, is a single chain of sites with the higher density 
of holes. That chain separates two AF domains with 
the opposite direction of sublattice magnetization. It is 
expected that the behavior of charge and spin excitations 
in the stripe of that kind resembles the physics of the t-J 
model (tJM) in one dimension (ID), the characteristic 
feature of which is the separation of spin and charge 121, 
Hoj l . We consider a different scenario according to which 
the stripe is formed when the DW which takes the form 
of a two-leg ladder is filled with holes. We also discuss 
the possibility of bond order formation in the ladder-like 
DW either on rungs or on legs. The isolated two-leg 
ladder is already an interesting object to analyze 0,^3. 
In copper-oxides, ladder-like arrangements of Cu atoms 
develop in the Cu02 planes containing line defects [T^ . 
Spin ladders display shor t range singlet correlations and 
gap in the spin spectrum 13, 14]. Systems with this type 
of behavior are often described by models of spin liquids 
or dimcr solids. 

As we have already mentioned, by analogy with the 
two-leg ladder we analyze in this paper the DW which 
consists of two chains that form the ladder. We also as- 
sume that spin singlets develop cither on legs or on rungs 
in the DW. The motion of a hole in the two-leg ladder 
is completely different in character than the motion of 
a hole in the AF background 0, Q, . The exchange 
interaction gives rise to formation of singlets on rungs. 
The hoping of the single hole between two sites along 
the same rung favors an even bonding state. The bond- 



ing state represents, by definition, a single fermion on 
the rung with a well-defined spin. The propagation of 
the hole along the ladder is related to a process during 
which the bonding fermion and the singlet on nearest 
rungs exchange their positions. Therefore, the motion 
of the hole in the ladder may be represented to lowest 
order by a tight binding model of spin-^ fermions 
The hopping parameter related to the propagation of the 
bonding fermionic state is reduced by a factor of 2 rel- 
ative to the value t of the bare hopping parameter. If 
the DW between two AF DWs has the form of a lad- 
der it is also natural to consider the formation of singlets 
on legs as a possible form of bond order. There exists 
some evidence that the AF order is formed in cuprates 
in hole-poor domains separated by hole-rich stripes 0]. 
It is believed that the phase of the staggered magnetiza- 
tion changes across each DW between AF domains by tt. 
The motion of a hole in the AF is governed by an energy 
scale ~ J which is lower than t or t/2. The hopping hole 
shifts spins and creates defects in the AF ordered state. 
The presence of defects brin gs a bout an increase of energy 
and confinement of the hole [ig . That confinement is not 
perfect because the defects may be annihilated by the ex- 
change term in the Hamiltonian which flips antiparallel 
spins on nearest neighbor (NN) sites. In our calculation 
we do not consider the coherent propagation of holes in 
the AF. This motion is mediated by the exchange term 
^ J < t, and it thus is a higher order process. This 
process may give rise to tunneling between stripes. We 
neglect effects related to it in the quantitative analysis, 
but we will later discuss consequences which it may have 
for the form of the stripe phase. 

Now, we temporarily concentrate on the anisotropic 
version of the tJM which is the tJ^ model (tJ^M) defined 
by the Hamiltonian H^, 

= -t ^(4.c,,. + H.c.) + jY, S:S^, (1) 

(id) (id) 

where {i,j) denotes a pair of NN sites. The Hilbert space 
is restricted to states representing configurations without 
doubly occupied sites. The tJM contains an additional 
term, 

H-H,, = J/2 ^(5+5- + SrS+), (2) 

which is needed in order to give the exchange interac- 
tion the isotropic form. We start with the discussion of 
hole motion inside the ladder-like DW with bond order 
on rungs. The physics of this motion is different in one 
aspect from the propagation of a hole in the ladder. The 
exchange interaction with the AF ordered domains acts 
on the DW as an effective staggered magnetic field. We 
may assume for convenience, that spins up occupy sites 
in domains on both sides of odd numbered rungs in the 
DW, as it has been presented in FigQIa), where numbers 
which label rungs are explicitly shown. An oval repre- 
sents a bond singlet. The bonding state on a rung is an 
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even combination 



1^(1) 



1 



m,a.L/ ' I m.a.Rl 



(3) 



of states l^'I^^Q l) l^m^a _r) representing a single 
fermion with spin a on the left site and on the right site in 
the m-th rung respectively, Fig^b) and (c). Apart from 
the m-th rung which is in the state i) or 1^'^^^^ 

all rungs are occupied by singlets. The hopping energy 
along the rung achieves its minimum value in the bond- 
ing state. Since the problem of the hole moving along 
the ladder-like stripe with bond order on rungs posses 
the axial symmetry, without any loss of generality we 
will later use in our considerations even combination of 
reflected states representing holes residing on both sides 
of the stripe axis. The bonding state Q is an example 
of such an even state. The meaning of the label "(1)" 
in Eq.Q will become clear later. Greek letters in FigQ] 
represent spins, a = —a and a circle denotes a hole (an 
empty site). For the matrix elements of the hopping term 
Ht in the < JM or the i J^M we have 



,(1) 



(4) 
(5) 



The contribution to energy from the magnetic part of 
the tJM from states and FigEIb) 

and (c), is higher if spins a and a are parallel and has a 
lower value if they are antiparallel. For parallel spins that 
contribution is higher by J/2 relative to the energy of the 
empty DW, Fig0Ja). For antiparallel spins a and a, the 
magnetic contribution to energy from the single fermion 
on a rung, Fig^b) or (c), does not increase relative to 
the energy of a rung occupied by the singlet. 

A state representing the single spin-up bonding 
fermion which propagates inside the ladder-like DW with 
the spin-Peierls order on rungs is given by the coherent 
sum, 




A2nk\\ 



(1) 



(/3o,T(fc||)l*^„^ 



-1,T' 



,T(fc||)l*ST))' 



(6) 



where fc|| is the momentum in the direction parallel to the 
stripe. It is probably worth to stress here, that the mo- 
tion of the hole is now restricted by definition to the DW, 
which is modeled by the two-leg ladder. Two different pa- 
rameters /3o,t(fc||) and /3e,t(fc||) have been introduced for 
odd and even numbers labelling rungs, because the en- 
ergy of the spin-spin interaction between the fermion in 
the DW and spins in the domains is different in the t J^M 
for fermions occupying even and odd rungs. L denotes 
the total length of the stripe. The matrix element of Ht 
between bonding states on NN rungs is Q, 



t/2. 



(7) 



By collecting all matrix elements we may deduce that 
the amplitudes /3e,f (fc||) and l3o.-\(k\\) obey the eigenstate 
equation, 



where 



H{k\\ 



-t + J/2 tcos(fc||)e 
tcos(fc||)e*''ii -t 



(9) 



The groundstate of H{k\\) is given by the following solu- 
tion of (jHJ 



/3o.T(fc||) = 



/3e.T(fc||) = 



^(|)2 + (iC0s(fc„)) 



2 _ I 
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-tcos(fc||)e*'=ii 



(10) 

(11) 



where denotes a real normalization constant. The 
groundstate eigenenergy e(fc||) is 



e(fc||) = -t + J/A ~ ■y/(J/4)2 + (tcos(fc||))2. (12) 

The interaction of the bonding fermion with the near- 
est spin in domains generates some effective magnetic 
field. Since this effective magnetic field which acts on 
the propagating bonding fermion is staggered, the effec- 
tive elementary cell gets doubled in the direction parallel 
to the DW. Thus, the first Brillouin zone (IBZ) of the 
stripe gets reduced by a factor of 2 and the quasiparticle 
band becomes split into a bonding band and an anti- 
bonding band. From now, e(A:||) in (jSJ will represent the 
energy of the bonding band which may be filled twice in 
the reduced BZ by propagating rung fermions with spin 
up or down. 

Now we proceed to discuss the motion of the hole in the 
direction transversal to the DW. The hole which enters 
an AF domain leaves a sequence of defects on its way. It 
must retreat to the stripe along the line formed by these 
defects, in order to repair them. The basis of states, 
which may be used to describe the hole motion along the 
DW consists of wavefunctions obtained during consecu- 
tive hopping in the domain, without retreats, of the hole 
which has started from a site in the DW ■ The creation 
of a next state belonging to this basis is mediated by the 
hopping term in the Hamiltonian. The recursion method 
is an approach useful to analyze a Hamiltonian H act- 
ing in a basis which consists of states created by applying 
consecutively that Hamiltonian to a given initial state |1) 
[l7| . This method has been already applied to discuss the 
hole motion in the Ising AF |l^ and to discuss the hole 

propagation along the narrow DW n m. Before we 

use that approach to the analysis of the composite hole- 
motion along the ladder-like DW we will shortly sketch 
basic principles of the recursion method. In general, any 
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Hamiltonian generates a basis of orthogonal states ac- 
cording to the following set of rules, 



|2> 



|n+l) 



(1|1> 
ln\H\n 



H\l) 

H\n) , 

{n-l\H\n) 
{n-\\n-iy 



n-l), 



(13) 



(14) 



where n > 2. An additional assumption is that none 
of states created according to the recipe (I14|l vanishes. 
The Hamiltonian is tridiagonal in the basis formed by 
states \n) and its eigenenergies may be found by analyz- 
ing the poles of the diagonal Green's function 511 (w) = 
1) which takes the form of a continued fraction, 



1 



511 (^) = 



(111) 



(i|g|i) 
(i|i> 



(2|2) 



(15) 



(i|i>('' 



<2|2> 



(3|3) 



As we have already mentioned, the hopping of holes in 
the ladder favors, in terms of energy, the even bonding 
states on rungs occupied by a single fermion, while the ex- 
change interaction favors the singlets on rungs occupied 
by two fermions. We choose a coherent sum of states 
representing the bonding state and the singlets on rungs 
© for the state |1) which we will apply in the recur- 
rence method. To use that state implies explicitly the 
assumption that the motion of the hole along the ladder- 
like stripe has two components. The motion inside the 
DW is analogous to the propagation between rungs of 
the bonding state, while the motion of the hole which 
has left the DW and has entered a domain is governed 
by processes related to the formation of spin polarons in 
the AF background 0, . By acting successively with 
the kinetic part of the t JzM on the state which repre- 
sents the hole created in the DW we create string-like 
states representing the hole that enters a domain from 
a given rung in the ladder-like DW. Here, by a string is 
meant a line consisting of defects created in the domain 
by the hopping hole which has started from a given site 
in the DW. A combination of states obtained in this way 
may be viewed as a spin polaron pinned to the ladder 
at the rung to which belongs the initial site from which 
the hole has commenced its motion. The recurrence ap- 
proach provides a convenient method to construct string 
states by applying the hopping term to the wave function 
©. Since ^ alone represents the coherent motion be- 
tween rungs of the bonding fermion, it is clear that such a 
construction takes into account the composite character 
of the hole motion along the stripe. 

The whole procedure looks in detail as follows. We 
use the state |^'^"'^''(fc||)) as the initial state |1) in the 

- I,T,(2) 



States |2) ee |^''^^^(fc||)), |3) = 



recurrence procedure 

,(3) 

■T . 

the hopping term and represent hopping of the hole in- 
side domains. As we have already stated, the initial state 



'^'X''{k\\))^ . . . are obtained by the consecutive action of 



|1) = I'J'j"'^'' (fc|| )) accounts for the component of the hole 
motion which is related to the propagation of the bonding 
rung-state inside the DW. The matrix element (1|H|1) in 
the continued fraction IjlSII can be deduced from the form 
of the Hamiltonian ((SJ 



,(1) 



X cos(2/cj|) 



H.c] 



2l/3o,T(fc||)l' 



(16) 



The Ising contribution to energy from the Hamiltonian 
Hz of i J^M in the previous expression is measured rela- 
tive to the energy of the empty DW. Since (3o,-\{k\\) and 
Pe.]{k\\) define the groundstate of JSJ, the matrix ele- 
ment (11611 reaches for that choice of /3oj(fc||) and (3e,^{k\\) 
its minimum value as the function of the parameters 
/3o,t(fc||) and /3e,t(fc||). This remark follows from the obvi- 
ous equivalence of jHll and (|16|l . We will later perform an 
additional variational analysis to demonstrate that the 
natural ansatz for |1) based on the solution of ||SJ) gives 
rise to qualitatively correct results. 

The hopping term moves the hole out of the DW deeply 
inside the domains. States obtained after the n-th hop 
of the hole contribute to the n-t-l-th state in the se- 
quence |1), |2), |3),... which is created when the re- 
currence method is applied to analyze string-like excita- 
tions in the ladder-type stripe. States depicted in Figs, 
md) and (e) contribute to the vector |2) = \^^^\k\\)) 
which appears in the definition of the Green's function 
(|15|l . The contribution to the Ising energy in the tJzM 
from a state obtained by one hop of the hole into the 
domains from a rung occupied by a single fermion is 
higher if the spin direction of the rung fermion and the 
direction of the nearest spin in this domains are par- 
allel, Fig^d). If these directions are antiparallel that 
contribution is lower because more bonds are occupied 
by antiparallel spins in the state created after a single 
hop of the hole, FigQIe). The hole which enters deeper 
into domains frustrates more AF bonds and its energy 
increases, Fig^f)-(h). After the first hop the hole has 
several choices between different directions of the next 
move inside the domain. There are exactly (z — 1) possi- 
ble directions of the next hop during which the hole does 
not return to the site occupied by it before the previous 
hop. z = 4 is the coordination number. Each next move 
may in general give rise to {z — 1) new states, examples of 
which have been depicted in Figsmf)-(h). The newly cre- 
ated states contribute to the next vector in the sequence 
|1) ^ |*«(A:||)), |2) ^ \^f{k\\)), |3) ^ |*f^(A:||)),... 
Scalar products which appear in (|15fl are given by, 



*«(fc,|)|vl/^^^(fc||)) 



(1)/ 



(17) 



(*^"'(fc||)|vl/("^(fc[|)) = (z-l)"-2t2(«-l)^ (18) 

where n > 2. The formula H17|l is merely the normal- 
ization condition for the amplitudes /3o,t and /3e.T- The 
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factor {z — in (|18|l represents the number of differ- 

ent paths along which the hole can move without retreats, 
while the factor t^Cn-i) appears there because the hop- 
ping term transforms the state \n) into the state |n-|- 1). 
The diagonal matrix elements of the Hamiltonian for the 
tJzM, between states = |n) representing holes, 

which have left the DW and have entered domains have 
their origin in the Ising term in that Hamiltonian and are 
given by, 

(*S")(fcii)ii7.,jvi/;")(fc||)) = 

X ^(" -^ + l/3o,T(fc||)p), 

(19) 

where n > 2. The formula H19|) has been obtained by 
counting the number of frustrated AF bonds. H19() de- 
pends on the proportion between the weights |/3o,t(fc||)P 
and |/3e.|(fc||)p because the number of frustrated bonds 
is higher if the hole starts from a rung occupied by a 
fermion with spin parallel to the nearest spin in domains 
and is lower if these spins are antiparallel. The energy 
of the propagating quasiparticle in the ladder-like stripe, 
ei(fc||) measured relative to the energy of the empty DW 
may be found within the recurrence procedure by looking 
for the lowest pole in (fT5|l . 

Strings formed by defects which have been created by 
the moving hole may be shortened in the isotropic sys- 
tem by the transverse part of the Heisenberg model which 
swaps antiparallel spins on NN sites. That process shifts 
effectively the hole or, to be more precise, the spin po- 
laron to a second or a third NN site. As we have al- 
ready mentioned, we mean by the spin polaron a state 
which is a combination of all string-like states created 
by the hopping hole which starts from a given site. We 
have also stressed that the propagating quasiparticle in 
the ladder-like stripe is the composition of the bonding 
fermion which moves inside the DW and the semicircular 
spin polaron. The latter characteristic manifests itself 
when the hole moves inside domains. Processes which 
give rise to the propagation of the polaron in the isotropic 
AF may bring about evaporation of holes from the single 
DW and it is likely that the single stripe is not stable 
in the doped isotropic AF described by the tJM. When 
stripes form an array, holes can not escape to infinity. 
In addition, it is clear that the system of stripes will 
be additionally stabilized by tunneling of holes between 
them. The full quantitative analysis of the stripe system 
and tunneling phenomena is beyond the scope of this pa- 
per. Here we concentrate on energy changes related to 
the internal structure of the single stripe. We will argue 
that the energy scale of these changes is higher than the 
energy scale of hole propagation in domains and tunnel- 
ing between stripes. Thus it is sufficient to consider the 
motion along the single stripe to draw some conclusions 
about the stability of the whole stripe system. The mo- 
tion of the hole inside the DW does not bring about any 
persistent changes in the spin structure, apart from some 



triplet excitations. It has been demonstrated that triplets 
play a minor role and may be neglected Q. Thus, the 
longitudinal component of the quasiparticle propagating 
along the stripe resembles a moving free fermion. That 
component which is related to the motion in the ladder 
of the bonding state on a rung, is crucial for lowering the 
energy of the stripe relative to the energy of the doped 
homogeneous AF. Since the energy scale for the polaron 
propagation in the AF is J, while the energy scale 
for the propagation in the DW (Fig^a)) of the bonding 
fermion is ~ t, in the lowest approximation we may ne- 
glect in the analysis of the stripe stability the processes 
which give rise to coherent propagation of polarons in the 
isotropic homogeneous AF. On the other hand, we con- 
sider the motion inside domains of holes which start from 
a given rung in the DW, retrace their steps and return 
to the initial site. That motion brings about a consid- 
erable change in the energy of the quasiparticle propa- 
gating in the stripe relative to the energy of the bonding 
fermion propagating in the ladder. This change is of or- 
der t. Thus we treat the propagating quasiparticle in 
the ladder-like stripe as a composition which consists of 
the bonding rung fermions and of the confined semicircu- 
lar spin polarons which are formed in AF domains near 
the DW. The reduction of the distance between stripes 
may bring about an energy change which is also of order 
^ t. Namely, the length of strings is effectively limited 
by the width of domains. We will take into account in 
our calculation the effect which the domain width has on 
the energy. 

We have already mentioned that the translational in- 
variance of the stripe in the direction parallel to its axis is 
broken by the AF ordering of domains. Thus the energy 
band formed by holes propagating in the stripe becomes 
split and the size of the ID Brillouin related to the mo- 
tion along the stripe is effectively reduced by a factor of 
2. The dotted fine in Fig. El represents the lower part of 
this band. It can be filled twice by quasiparticles which 
are compositions of the spin polaron and of the bond- 
ing spin-up or spin-down fermion on a rung. A stripe 
obtained by filling the ladder-like DW with holes is dif- 
ferent in this respect from the stripe obtained by filling 
with holes the narrow DW 0. The quasiparticle band 
which forms for the latter type of stripe is represented 
by the dash-dotted line in Fig|21 It can be filled only 
once because a spin-less holon is now the component of 
the quasiparticle propagating along the stripe 0. The 
dispersion of that band has been also found by means of 
the recurrence method. The energy is measured in both 
cases relative to the energy of the empty DW. Irrespec- 
tive that the bottom of the energy band for the stripe 
which is formed when the ladder-type DW is filled, lies 
higher than the bottom of the energy band for the stripe 
which is formed when the narrow DW is filled, filling 
the ladder-type DW with holes is a more effective way of 
lowering the energy for higher levels of doping because 
the quasiparticle band is flatter in this case and it can 
be filled twice. The dashed, doubly-dotted, straight line 
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represents the energy of a hole confined in the anisotropic 
AF described by the tJ^M. 

Until now we have assumed that the amplitudes 
/3o,t(^||) ^iid /3e,T(fc||) which define the state |1) = 

are solutions of the eigenstate equation ((HJ. 
It might be interesting to analyze if that natural choice, 
which is optimal in terms of lowering the matrix element 
(|16|l that describes motion of the hole inside the DW, 
provides a good approximation to the wave function of 
the hole propagating along the stripe. We will check how 
good is the initial choice of |1) — \'^''^^\k\\)) by varying 
the amplitudes /3o,T(fc||) and (3e,^{k\\) and looking for the 
minimum of the lowest pole in H15|l considered as a func- 
tion of these parameters under an additional condition 
that they are normalized according to the relation 

l/3o,T(fc||)l' + l/?e,T(fc||)P = l. (20) 

Now, |/3o,t(^||)P (|/3e,T(fc||)P) represents the weights with 
which polarons pinned to odd (even) numbered rungs 
contribute to the quasiparticle. The weight l/^o,!^^!!)^ 
may be different from |/3ej(fc||)p because the contribu- 
tion to the Ising part of the energy in the tJ^M is lower 
if the hole enters the domain from a bond occupied by a 
fermion with spin antiparallcl to the nearest spin in the 
domains. We have initially used fixed values of (3o,-\{k\\) 
and /3e,t(fc||) in agreement with the original scheme of 
the recurrence method in which the form of the initial 
state is fixed. Now we combine the recurrence method 
with the variational approach in order to check the cor- 
rectness of the initial ansatz for |1) and the sensitiv- 
ity of results to the change of the ansatz parameters 
Po,i{k\\) and /3e,-f (fc||). The dashed line in Fig|21represents 
ei{k\i) obtained by optimizing the amplitudes Po,iik\\) 
and /3e,t(fc||). Despite some differences between the dot- 
ted and dashed curves we may conclude that the initial 
ansatz for |1) = |\I'j"'^-'(fc||)) gives rise to the qualita- 
tively correct description of the energy dispersion e;(fc||) 
of the quasiparticle propagating in the ladder-like DW. 
The optimization of parameters Po.i{k\\) and /3e.|(fc||) 
lowers the energy of the propagating quasiparticle, but 
the change is not dramatic which indicates that even 
the original ansatz has been reasonable and the whole 
method is not very much sensitive to the change of de- 
tails. 

It is worth to point out, that the combination of the 
recurrence approach with the optimization of the pa- 
rameters /3o,t(fc||) and /3e.|(fc||) is still an approximation. 
States contributing to a given vector \n) and related to 
strings of length n — 1, obtained by hopping of the hole 
which starts from even and odd rungs are created in the 
same step in the recurrence procedure. Thus, the ra- 
tio between weights of strings with different lengths is 
the same for both kinds of strings independently of their 
length, which is not necessary optimal in terms of energy 
lowering. Since the difference between Ising energies of 
string states obtained by hopping of the hole which has 
started from an even or an odd rung does not change 



much with the length of strings that approximation defi- 
nitely should not influence the final result in a qualitative 
way and thus it is acceptable. 



III. HOLE MOTION ALONG THE 
LADDER-LIKE STRIPE WITH SPIN-PEIERLS 
ORDER ON LEGS 

In this part of the paper we consider a different under- 
lying spin structure of the bond-centered stripe. It is not 
clear what kind of bond order is energetically favorable 
when a ladder-like DW is filled with holes. An obvious 
optional magnetic structure is the spin-Peierls order on 
legs. Figl^Ja) depicts that structure. We sketch now the 
calculation of the energy dispersion for the quasiparticle 
moving in such a spin background. The hopping of a 
hole along a rung gives rise to an exchange of positions 
between a single fermion and a singlet on parallel legs. 
The matrix element of the hopping term in the Hamilto- 
nian which couples the states depicted by FiglSfb) and 
(c) or couples the states depicted by FiglSJ^d) and (e) is 
t/2. Since that matrix element is positive, the lowest en- 
ergy state with a given momentum along the stripe may 
be constructed by summing odd combinations of states, 
which are images obtained by reflection in the axis of 
the stripe. Two examples of reflected states are shown 
in Fig|2Ib),(c), and in Fig|2Id),(e). In a similar way as 
for stripes with the spin-Peierls order on rungs we define 
l^'n^u.ct) as the normalized odd combination of states de- 
picted by FigEfb) and by FigE^c). denotes the 
odd combination of states depicted by Fig|3|[d) and by 
FigEIe). A wave function which we use as the starting 
point in the recurrence procedure applied to the prob- 
lem of the stripe with spin-Peierls order on legs may be 
defined as 

+ 0uAh)K%))- (21) 

The parameters Pi^^ and (3u.a- obey the normalization 
constraint 

+ = 1. (22) 

In an analogous way as before, we may deduce that the 
diagonal matrix element for that initial state is 

{¥^\k\^)\HuA^^^Hk\\)) = -l-mAh)(3iAk\\) 

+ /?,%(/c||)/3„,.(/c||)] 

+ ^IPuAh^- (23) 
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By minimizing the matrix element 1)23(1 with respect to 
the parameters (3i^„, and /3„_o- "we may determine their 
values which can be later used as an ansatz in the re- 
currence procedure or alternatively we may look for a 
better approximation by performing the minimalization, 
after the lowest pole of 1)21(1 has been found. Re- 

lations presented below may be obtained by carrying out 
the recurrence analysis in a way similar to that which 
was applied to the problem of the ladder-like stripe with 
the spin-Peierls order on rungs 



(Vl/W(fc||)|vl/W(fc||)) = 1, 
(vl.(«)(fc||)|*(«)(fc||)) = 



1-2, 2{r. 



(24) 
(25) 



and 



X (n-l)J, (26) 

where n > 2. By looking for the lowest pole of gii{uj), 
Eg. ((15(1 . and optimizing that parameter with respect to 
I3i,a-{k\\) and /3u,cr(fc||) we may find the energy dispersion 
ei(A;||) of a hole- like quasiparticle propagating along the 
ladder-like stripe with the spin-Peierls order on legs. The 
solid line in Fig|5] represents that dispersion. As before, 
ei(fc||) is measured relative to the energy of the undoped 
DW forming the stripe upon hole filling. The solid line 
in Fig[21 actually simultaneously represents the disper- 
sion obtained with the parameters Pu,crik\\) and /3;.g-(fc|| ) 
derived by minimalizing only the matrix element ((23(1 . 
The lack of difference between two approaches may be 
attributed to the fact that the matrix element of the 
Hamiltonian for higher states created in the recurrence 
procedure does not depend on values of the parameters 
/3u,o- (fc||)) and /3i^cr(fc||)- Thus we may conclude that in 
the case of the stripe with the spin-Peierls order on legs 
already the ansatz wave function \"^lj\k\\)) derived by 
minimizing the matrix element 1(23(1 gives rise to the sat- 
isfactory description of the hole propagation along that 
stripe. 



IV. BOND-ORDERED ANTI-PHASE DOMAIN 
WALL IN THE ANTIFERROMAGNET 

In this section, as all over this paper, we again view 
the stripe as a ladder-like DW, to some extent filled with 
holes. Now, we discuss the undoped system and ana- 
lyze the energy of the DW with a given internal mag- 
netic structure by means of the second-order perturba- 
tion theory applied to the Heisenberg model. First we 
discuss the DW with the dominating spin-Peierls order 
on rungs. Our starting point is an unperturbed Hamil- 
tonian Hq with a groundstate l^o), shown in Fig^Ja) 
which depicts two AF domains separated by a column 
formed by singlets. 



E 

(m,n) 



(27) 



where refers to all pairs of NN sites belonging to 

domains and (m, n) refers to all pairs of NN sites belong- 
ing to rungs in the ladder-like DW. As before, ovals in 



FigEIa) denote singlets on rungs. The operator s\j 



V2 



V t t 



(28) 



creates a singlet on an empty pair of sites i, j- It is 
worth to stress there that the whole exchange interaction 
of sites which belong to domains with sites belonging to 
ladder-like DWs is treated here as the perturbation. In 
the second-order perturbation theory the energy of an 
eigenstate |$) of iJ is, 

y^ ($(")|gi|^f))(vl/f)|gi|ci>(")) ^ ^^g^ 
Ei — Eq 

For the problem of the DW in the two-dimensional AF, 
Hi = H — Hq defines the perturbation term in the Hamil- 
tonian H of the Heisenberg model on the square lattice, 
1'!'^°^) are excited states of Hq with eigenenergies Ei and 
Eq is the groundstate energy of Hq. |<i>(°)) is its ground- 
state. Figures0Jb)-(g) depict typical examples of excited 

states l^'l"'') which contribute to 1(29(1 . A zig-zag line rep- 
resents a "broken" bond with a higher contribution to the 
energy of the Ising model relative to the energy of the 
Neel state. A letter o= x,y, z inside a rectangle denotes 
a triplet state created by the operator 



V2 



(30) 



acting on a an empty pair of NN sites. 

The difference between the energy of the ladder-like 
DW of length L and the energy of the homogeneous 
ground state of the quantum AF in 2D (Heisenberg 
model) which can be calculated in the second pertur- 
bation theory by collecting contributions to ((29(1 from all 
excited states is. 



5E, 



13 

40^ 



-JL. 



(31) 



The groundstate energy of the homogeneous AF in 2D 
which has been used to derive the value of SEi, in l(31() 
has been also obtained in the second-order perturbation 
theory. The Ising model and the Neel state have been ap- 
plied in the derivation as the unperturbed Hamiltonian 
and its groundstate, respectively. We would like to point 
out that in order to analyze eigenstates of the Heisenberg 
model representing the homogeneous state and different 
kinds of DWs we must choose different parts of Hamilto- 
nian as unperturbed terms. 

We proceed now to discuss the energy of the narrow 
DW between two AF domains with the opposite stag- 
gered magnetization. The domains are in this case in 
direct contact with each other, as shown in Fig0Jh). An 
increase of the energy due to the formation of the narrow 
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DW is lower than for the wider DW depicted by FigQfa), 



6 



JL. 



(32) 



The quantity 5En has been obtained also by means of the 
second-order perturbation theory applied to the Heisen- 
berg model. Quantum fluctuations in the isotropic AF 
reduce the increase of the energy caused by the DW. 
We may draw this conclusion by noticing that in the ex- 
tremely anisotropic case of the Ising model the energy of 
the ladder- like DW, Fig0Ja), is JL while the energy of 
the narrow DW, Figl^h), is \JL. 

We proceed now to discuss the energy of the wide DW 
with the dominating spin-Peierls order on legs. The un- 
perturbed Hamiltonian may be written in this case also 
in the form given by H27|) . but with pairs (rn,n) refer- 
ring to legs on which singlets are formed, as it has been 
presented in Fig0Ji). Again, by analyzing excitations 
created in the groundstate of Hq by the remaining part 
of the Hamiltonian we may get an approximation to the 
energy of the DW with the spin-Peierls order on legs. It 
turns out that this energy is j^JL which is lower than 
for the DW with the spin-Peierls order on rungs. This 
lowering may be mostly attributed to the existence of a 
bigger number of accessible excitations than in the pre- 
vious case. An example of a new kind of excitation has 
been presented in FigQJj). The energy of the DW with 
the spin-Peierls order on legs is JL in the Ising limit. This 
value is the same as for the DW with the spin-Peierls or- 
der on rungs. Again we notice that quantum fluctuations 
diminish the energy difference between the wide and the 
narrow DWs. We will show that this reduction is crucial 
for the stability of the ladder-like stripe. 



V. STABILITY ANALYSIS 

ei , the energy per one doped hole of the stripe obtained 
by hole filling the ladder-like DW with the spin-Peierls 
order on rungs is given by the obvious formula. 



6/ 



-ft dkn 

TT <J —kw II 



5Ei 



(33) 



where 5Ei is the energy of the DW per one lattice spacing 
along the stripe, e/ and 5Ei are measured relative to the 
energy of the undoped AF. The linear filling of the stripe 
is 



1 



kp- 



dkn 



2kF 

TT 



(34) 



(|33|l and H34II are true for the doping level lower than 
1 hole per site, only. Formulas which may be applied 
to derive the energy per one doped hole and the linear 
filling parameter for the stripe with the spin-Peierls order 



on legs are similar. A main difference is that the energy 
band has now its minimum at the point 7r/2, 



ei 



and 



2-fei 



ez(fc||)dfc|| +5Ei 



TT J —TV 



-Tr/2+kF 



/2 



dku 



1 

TT JTT/2 — kp 



dku 



(35) 



<5 =-( 



-TT/2 + kp 



dku 



-tt/2 



/ dfc||) = f^. (36) 

JTT/2-kF 



The site-centered stripe may be formed by filling with 
holes only one of two vertical chains in the middle of 
FigEfh). The holon band which is formed for the site- 
centered stripe obtained by filling the narrow DW in 
FigEIh) can be be filled only once because holons are 
spin-less objects 0. Thus, the energy per one doped 
hole e„ is given by 



1 

2tT J TT 



TT + kp 

kp 



i(fc||)(ifc|i -I- SE„ 



h f dk 

In J -n — kp 



(37) 



where 5En is the energy of the narrow DW per one lattice 
spacing along the stripe. The linear filling parameter 5n 



5l 



1 

2^ 



TT + kf 



dku ~ 



-kt 



(38) 



Fig|Sl depicts the total energy per one doped hole in the 
ladder-like DW (dotted line, dashed line and solid line) 
and the total energy per one hole in the narrow DW 
(dash-dotted line) as a function of linear filling jj'^^-j in 
the extremely anisotropic limit of the tJ^M. Dash-doubly 
dotted horizontal line represents the energy of the sin- 
gle spin polaron confined in the homogeneous Neel state. 
The dotted line represents the energy e/ for the ladder- 
like stripe with the spin-Pcicrls order on rungs derived 
from the formula (|33|l in which we have applied the quasi- 
particle energy ez(A:||) estimated by means of the calcula- 
tion in which the parameters /3e,t (fc|| ) and /3o.t(^||) have 
been determined by the ansatz that gives rise to 
the minimal value of the the matrix element H16(l . The 
dashed line has been obtained for the same type of stripe 
by minimizing the value of the lowest pole in (|15|) as a 
function of these parameters. The solid line refers to the 
wide stripe with the spin-Peierls order on legs. The en- 
ergy is measured in all cases relative to the energy of the 
Neel state. Motivated by estimates based on experimen- 
tal data we have chosen J/t = 0.4. The stripe obtained 
by filling the narrow DW becomes stable slightly above 
quarter-filling in the limit of extremely anisotropic ex- 
change interaction of the tJ^M. The energy per one doped 
hole of the stripe obtained by hole filling the ladder-like 
stripe with the spin-Peierls order on rungs is lower than 
the energy of the hole confined in the Neel AF for the 
linear filling above Sf ~ 0.8. The energy per one doped 
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hole for a ladder-like stripe with the spin-Peierls order on 
legs is higher than the energy of the hole confined in the 
anisotropic AF for the whole range of the linear filling pa- 
rameter. Since the energy of the ladder-like DW, is much 
higher than the energy of the narrow DW, the ladder-like 
stripe has higher energy than the narrow stripe up to the 
linear filling parameter sf ~ 1.0. Nevertheless, for val- 
ues of Sf, . < 0.5, the hole-filling of the ladder-like DW 

l{n) ~ ' ^ 

is a more effective way of lowering the energy relative 
to the energy of the empty DW, than the hole-filling of 
the narrow DW, because the quasiparticle band for the 
ladder-like stripe is flatter, lies below the average energy 
of the quasiparticle band for the narrow stripe and it can 
be filled twice. Fig|H| depicts the difference between the 
energy per one doped hole in the ladder-like stripe and 
in the narrow stripe as a function of doping <5|'j^j. The 
dotted (dashed) line again represents results of the cal- 
culation for the ladder-like stripe with the spin-Peierls 
order on rungs with non-optimized (optimized) parame- 
ters /3e,^{k\\) and /3o,T(fc||)- The solid line refers to results 
of a calculation performed for the ladder-like DW with 
the spin-Peierls order on legs. We stress that energies 
are measured now relative to energies of the empty DWs 
because contributions to ei and e„ from SEi and SEn 
have been omitted. Thus, this figure solely demonstrates 
the effect which on the energy has the filling of bands 
depicted in Fig|21 The hole filling of the ladder-like DW 
with the spin-Peierls order on legs seems to be a less effec- 
tive way of lowering the energy than the hole filling of the 
ladder-like DW with the spin-Peierls order on rungs. By 
looking at formulas H33|l - I|38(l we notice that the reduction 
of the difference between DW energies of the narrow DW, 
SE^ ^ and the ladder-like DW, 6Ei in §^ and ^ 
may give rise to the stability of the ladder-like stripe at 
finite doping. As we have already shown that difference is 
reduced in the isotropic AF and one may expect that the 
ladder-like stripe may be more favorable in the system 
described by the tJM. It turns out, that for the linear 
filling above (jj'^^^^ ^ 0.65 the ladder-like stripe becomes 
stable in the isotropic system, Fig[7| In the calculation 
performed to derive data which have been used to plot 
FigEI we have applied in formulas , and the 
DW energy per site, 6Ei, and SEn obtained for the fully 
isotropic AF. The ladder-like DW with the spin-Peierls 
order on legs has lower energy at a certain range of linear 
doping above that value. At higher values of the doping 
parameter, the DW with the spin-Peierls order on rungs 
is more stable. The stability range of the spin-Peierls or- 
der on legs in the DW slightly depends on details of the 
calculation. We have noticed this by comparing curves 
obtained by using the optimized and non optimized pa- 
rameters /3e,t(fc||) and /3oj(fc||). Notwithstanding that 
our results are somewhat sensitive to the applied calcula- 
tion method, it seems that the bond-centered stripes with 
the spin-Peierls order become more stable than the site- 
centered stripes for the linear filling parameter slightly 
above half filling. Despite that it is not possible to find 



within our approach the exact value of the doping pa- 
rameter at which the bond-centered stripes become more 
stable than site-centered stripes, it clear that spin fluctu- 
ations mediated by the transversal part of the Heisenberg 
model stabilize a system of stripes which are formed by 
hole- filling of ladder- like DWs with spin singlets on rungs 
or legs. On the other hand, it seems that a single stripe 
of that kind is not stable because holes which propagate 
in the isotropic AF will easily evaporate from a single 
DW. A system of such stripes is stabilized by tunneling 
of holes or hole pairs between stripes. Detailed analy- 
sis of such phenomena which seem to change the total 
energy to much lesser extend than hole kinetics along 
the stripe is beyond the scope of this paper. We ex- 
pect that the energy change which is brought about by 
the propagation of the hole through domains does not 
influence in a decisive way the shape of a stable stripe 
system, because hole hopping through the homogeneous 
AF background involves the exchange term which is pro- 
portional to J < t. There seems to exist in a system 
of stripes a mechanism which may bring about a change 
of energy which is related to hopping and does not in- 
volve directly the exchange term. In a system of stripes 
the length of strings created by holes which hop in do- 
mains is to some extend limited by the width of these 
domains. This limitation influences the band renormal- 
ization which is brought about by intrusions into domains 
of the hole moving along the DW. Neutron scattering ex- 
periments suggest that the width of AF domains is 2-3 
lattice spacings for the doping level 1/8. Now, we are go- 
ing to check how the energety of stripes is changed by the 
stripe-length limitation induced by the finite spacing be- 
tween stripes. The same quantities as in FigElhave been 
depicted in FiglHl The only difference is that to derive 
energies presented in Fig|Sl we have limited the length 
of strings to 3 lattice spacings. By the string length is 
meant here the number of hops needed to create a given 
string state from a state representing the single hole on a 
bond. By analyzing the construction of the sequence of 
states |1), 1 2)... by means of the recursion method applied 
to the stripe problem we easily infer that the state \n) 
represents strings of length \n — 1). Thus, by stopping 
the recurrence procedure at the n + 1-th step we may 
impose the restriction of the string length to n lattice 
spacings. We notice that the effective radius reduction 
which takes place for semicircular quasilocalized polarons 
which are formed in domains brings about a shift of the 
lower limiting value of the doping parameter, for which 
the bond-centered stripe is stable, towards a lower level of 
linear filling. Energy values depicted in Fig|51 have been 
obtained under a more restrictive assumption that the 
maximal value of string length is 2 lattice spacing. We 
conclude that the narrowing of domains is favorable for 
bond-centered, ladder-like stripes with the spin-Peierls 
order on legs. 

We do not assume that nearest DWs are impenetra- 
ble for a hole which has been created at a given DW. 
The hole may penetrate a nearest DW without rising the 
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exchange energy only if at a certain stage of this pro- 
cess the transversal part of the Heisenberg Hamiltonian 
is applied to an intermediate state. The transversal term 
is proportional to J. We believe that omitting such a 
process does not change final conclusions which we draw 
from this calculation. 



VI. DISCUSSION AND OUTLOOK 

The crucial question which needs to be answered, be- 
fore the puzzle of high temperature superconductivity 
and a non-Fermi liquid behavior observed in cuprates 
is solved, is how the Mott insulating antiferromagnet is 
transformed, upon doping, into the unconventional metal 
and the superconductor. It is clear that properties of the 
doped AF are determined by the competition between 
the kinetic energy, which favors formation of a metallic 
state with well defined itinerant quasiparticles, and the 
exchange interaction which favors magnetic or charge- 
ordered states that are insulators or poor conductors. 
A fully phase-separated state with one magnetically or- 
dered hole-poor region is favorable in terms of lowering 
the exchange energy but not in terms of lowering the 
kinetic energy. A compromise between two seemingly 
opposite tendencies to localize moments and to form a 
broad quasiparticle band is reached by means of micro- 
phase separation which takes the form of the stripe phase. 
Much effort has been made to develop theories of stripe 
states in doped AFs |^ 0, 120I l2l| . The experimental de- 
tection of stripes is far from being trivial. Their shape 
has not been identified yet with full certainty. Experi- 
mental evidence of micro-phase separation and its the- 
oretical interpretation is not always unique The 
width and the magnetic structure of stripes, defined as 
hole-rich regions which play the role of DWs for magnet- 
ically ordered regions, belongs to issues that are still un- 
der debate. The concept of magnetically ordered states 
in doped AFs, have been recently confronted with the 
idea of bond-ordered states. The classification presented 
in an article (ij suggest that magnetic and paramagnetic 
orders are optional types of ordering, which do not co- 
exist. A similar confrontation was already performed in 
an earlier paper in which energy changes upon dop- 
ing in the columnar spin dimerized phase and the Neel 
state have been compared. In some other articles 
the quantum critical point between the Neel state and 
the columnar dimerized phase has been examined. It is 
worth to emphasize that in all these publications it was 
assumed that dimers in the columnar phase were cover- 
ing the whole square lattice. We have presented in our 
paper a different point of view, according to which the 
coexistence of bond order with AF order is possible due 
to micro-phase separation. The doped antiferromagnet 
favors the magnetic order with nonvanishing value of the 
sublattice magnetization in hole-poor parts of the sys- 
tem, while the bond order develops in hole-rich regions 
which take the form of unidirectional stripes separating 



AF domains. The formation of the spin-Peierls order in 
stripes which play the role of DWs is favorable, because 
holes may propagate more freely in such a spin back- 
ground and the kinetic energy is substantially lowered in 
this way. Within the scenario of coexisting bond, AF 
and charge orders we have considered a single DW be- 
tween two AF domains with opposite directions of the 
sublattice magnetization. This DW takes the form of the 
two-leg ladder with the spin-Peierls order either on legs 
or on rungs. We have assumed that the motion of holes 
inside the DW is governed by a process during which a 
singlet on a bond and a single fermion on a neighbor 
bond exchange their positions. This location-exchange 
process is of first order because it is mediated by a single 
action of the hopping term. Thus it lowers effectively the 
kinetic energy. We have also assumed that a strong ten- 
dency towards confinement manifests in the motion of a 
hole which has left a stripe (a DW) and has entered one 
of two AF domains which are separated by the stripe. 
This quasiconfinement is related to the formation of so 
called strings or paths along which the hole has been 
hopping to a given site. The presence of strings which 
consist of defects in the AF structure and which are cre- 
ated by the moving hole, gives rise to the energy increase. 
The tendency towards confinement may be attributed to 
this energy rise. Holes, by following exactly the same 
path which they have taken to desert the stripe may re- 
turn to it without leaving any defects in the magnetic 
structure of domains. The motion of the quasiparticle 
along the stripe may be viewed as position exchange of 
a bond fermion and of a bond singlet between nearest 
bonds. That motion is strongly renormalized by hole in- 
trusions into AF domains. After calculating the band 
energy of such a propagating complex quasiparticle we 
have deduced that the system of bond-centered, ladder- 
like stripes is energetically more stable than the system 
of site-centered stripes in the isotropic AF at the linear 
hole-filling slightly below 0.5. We have also taken into 
account in our analysis that the form of bond order in 
the ladder-like stripe, treated as a DW, strongly influ- 
ences the value of the exchange energy. We have also 
argued that processes related to the coherent hole prop- 
agation in the AF background and to the hole tunneling 
between stripes are of lower order and that they change 
the energy to lesser extent. There exists some experimen- 
tal evidence, that the linear filling of stripes observed in 
cuprates is about 0.5. Since, in the results of our calcu- 
lation, the limiting value of the linear filling parameter 
for which the bond centered stripe becomes more stable 
than the site centered stripe is not far away from that 
number, we may draw a general conclusion from our cal- 
culation, that the competition of different order types 
may inffuence properties of moderately doped cuprates. 
The outcome of our analysis also demonstrates that the 
system prefers formation of singlets on bonds parallel to 
stripe axes. Such direction of bond order in stripes has 
also been observed in results of exact diagonalizations 
[25I |2^ . From the comparison of Figs[7| |H1 and M we 
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may draw the conclusion that bond centered stripes are 
formed for higher levels of global doping at which the dis- 
tance between stripes is shorter. The narrowing of AF 
domains brings about the shortening of the average string 
length and the increase of weight for states representing 
holes inside DWs. It seems that the motion of holes in- 
side the bond centered stripe with the spin-Peierls order 
is more effective in terms of lowering the energy than the 
motion of holes inside the site-centered stripe. The influ- 
ence of this factor becomes more pronounced if distances 
between stripes are small. We have observed that quan- 
tum fluctuation reduce the energy difference between the 
bond-centered DW and the site-centered DW in the 2D 
AF. Since the limiting value of the linear filling parame- 
ter at which the bond-centered stripe becomes stable lies 
in FigO below the number 0.5 which refers to the global 
doping level about 1/8 and the width of AF domains 
about 2-3 lattice spacings, we may deduce that the spin- 
Peierls order is very likely for this and higher amount of 
doping. 

Bond order has been already discussed in the context of 
stripes in some earlier papers. In one of these articles 
it has been demonstrated that the spin columnar dimer- 
ized phase is stable at a certain range of global doping. 
The author of this paper has suggested that such a state 
develops in cuprates as the stripe phase. It seems that 
it is desirable to consider some additional modulations 
of the spin structure than dimer columns in order to ac- 
count for the doping-dependent spacing between stripes. 
Magnetic excitations in the stripe phase have also been 
discussed in the literature. Spin wave theory of these 
excitations, based on a scenario of site-centered stripes 
separating AF domains 27] does not capture the high en- 
ergy part of the spectrum measured by means of neutron 
scattering |j3| . Assuming spin dimerization and formation 
of bond centered stripes authors of two recent preprints 
have been able to formulate a theory which correctly de- 
scribes the measured spectrum 28]. On the other hand, 
they have assumed that hole-rich regions taking the form 
of two- leg ladders do not participate in spin dynamics, 
which is hard to justify because at the linear filling 1/2, 
3 sites per 4 are occupied by spins in such structures. 

According to our scenario the bond order dominates 
in hole-rich stripes which play the role of DWs for hole- 
poor domains with more pronounced AF correlations. It 
is energetically favorable for holes to fill stripes because 
there is no tendency to confine hole motion in two-leg 
ladder-like structures with the spin-Peierls order. AF 
correlations dominate in hole-poor regions which are ef- 



fectively magnetically coupled via rungs occupied by a 
single fermion. If the spin direction of this fermion is 
antiparallel to the direction of spins lying in domains 
on both sides of the rung, by hopping forth and back 
along the rung, the fermion is able to lower the contribu- 
tion to the exchange energy on both sides of the DW. A 
similar mechanism is known to operate in site-centered 
stripes 29] . The two-leg ladder in the effective staggered 
magnetic field was already discussed, in the context of 
stripes, in some previous pubhcations 113, but authors 
of these papers either did not consider bond order inside 
the ladder or did not consider the formation of strings in 
domains and hole intrusions into them. 

The experimental identification of bond order in 
stripes and its shape is a difficult task. There is no di- 
rect method to achieve that goal. Important information 
is provided by inelastic neutron scattering experiments 
p]. On the other hand results of theoretical analyzes 
|28j ] seem to indicate that the shape of measured spectra 
is not very much sensitive to the form of the underlying 
bond order. Thus, it seems that by means of neutron 
scattering measurements alone is it not possible to iden- 
tify the exact form of the bond order which is formed in 
stripes. The interpretation of angle resolved photoemis- 
sion experiments may be additionally helpful in this task, 
because the formation of a bond ordered state influences 
the electronic structure of any system The single- 

particle spectral weight in the stripe phase was discussed 
in several theoretical papers plLISll ]. On the other hand, 
the impact which the shape of hypothetical bond order 
in stripes, understood as hole-rich DWs between AF do- 
mains, may have on spectral properties of cuprates has 
not been yet thoroughfuUy analyzed. Results of such an 
analysis performed by us confirm conclusions which we 
draw in this paper and will be reported elsewhere. 

Much evidence has been accumulated that in some 
cases the charge and spin inhomogeneity observed in 
cuprates has purely 2D character [33] ■ Such modulations 
seem to be also consistent with the scenario of fluctua- 
tions between some bond ordered states [3^ . An analysis 
of this category of states is beyond the scope of the dis- 
cussion presented in our paper. 
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FIG. 1: Graphical representation of some string states created 
by the hopping inside AF domains of a hole created at the 
ladder-like DW with the spin-Peierls order on rungs, a is 
equivalent to —a. 



2 . 6 



« 



-2 . 8 Li 








7r/2 



TT 



k 



FIG. 2: Energy of a single hole confined in the homogeneous 
Ising AF (dash-doubly dotted line), energy of a quasiparticle 
propagating in the ladder-type stripe with singlets on rungs 
(dotted line and dashed line), energy of a quasiparticle prop- 
agating in the ladder-type stripe with singlets on legs (solid 
line) and energy of a quasiparticle in the narrow stripe(dash- 
dotted line). Energy is measured relative to energy of the 
homogeneous Ising AF or an empty DWs. fcy is the momen- 
tum along the stripe. Only the symmetric halves of bands for 
positive fc|| have been presented. For the spin-Peierls order 
on rungs in the ladder-like stripe, hero and below, the dotted 
(dashed) line refers to a calculation with non-optimized (opti- 
mized) amplitudes /3e,T(fc||) and Po,t{k\\) of states representing 
transversal spin polarons attached to even and odd numbered 
rungs. J/t = 0.4 as everywhere. 
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FIG. 3: Graphical representation of some string states created 
by the hopping inside AF domains of a hole created at the 
ladder-like DW with spin-Peierls order on legs. 
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FIG. 4: Graphical representation of the DW which takes the 
form of a two- leg ladder with bond order on rungs, (a); and 
some excitations in this DW, (b)-(g). Graphical represen- 
tation of the narrow DW representing two directly merged 
anti-phase domains which are not separated by an intermedi- 
ate region with additional order, (h). Magnetic structure of 
a stripe formed by hole-filling of the DW with bond-order on 
legs, (i), and an excitation created in that structure by the 
transversal part of the Heisenberg model, (j). 
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FIG. 5: Cj, the energy per doped hole relative the energy of 
the homogeneous Ising AF for the ladder-like stripe with the 
spin-Peierls order on rungs, dotted lino and dashed lino, and 
for the ladder-like stripe with the spin-Peierls order on legs, 
solid line. e„, the energy per doped hole for the site-centered 
stripe, dash-dotted line, as a function of linear doping jj'^^j . 
The dash-doubly dotted line represents the energy of a hole 
confined in the homogeneous system. 



FIG. 6: Difference between the energy per doped hole in the 
laddcr-liko stripe and in the narrow stripe as a function of 
doping parameter Sf^^^^ . Since contributions to e; and e„ from 
6Ei and SEn have been neglected, energy is measured here 
relative to energies of empty DWs. The dashed (dotted) line 
represents that difi^erence calculated for the DW with the spin- 
Pcicrls order on rungs and the (non) optimized parameters 
/3e,T(fc|l)i and /3o,T(fe|i), while the solid line refers to the DW 
with the spin-Peierls order on legs. 
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FIG. 7: Energy per doped hole relative to the energy of the 
homogeneous isotropic AF for the ladder-like stripe, ei, dotted 
line, dashed line, and solid line and for the narrow stripe, e™, 
dash-dotted line, as a function of linear doping ■ 
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FIG. 8: The same quantities as in Fig|7| calculated under the 
assumption that the length of strings is restricted to 3 lattice 
spacings. 
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FIG. 9: The same quantities as in Fig|3 and |H1 calculated 
under the assumption that the length of strings is restricted 
to 2 lattice spacings. 



